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, Abstract. The components of the wide binary HD 219542 were recently found to differ in metallicity by about 0. 1 dex (Gratton 

CO . et al. l2001bl . In this paper, we present the results of 2 years of high precision radial velocity monitoring of these stars performed 

' at the Telecopio Nazionale Galileo (TNG) using the high resolution spectrograph SARG. No indication for radial velocity 

variations above the measurement errors (~ 5 m/s) was found for the metal richer component A. This allows us to place upper 
^ l' mass-limits for planets around this star. HD 219542 B instead shows a low amplitude variation with a 112 day period at a 

confidence level of ~ 96 - 97%. This might suggest the presence of a Saturn-mass planet, although it is still possible that these 
variations are due to moderate activity of the star. Tests based on variations of bisectors, stellar magnitude and line equivalent 
widths were inconclusive so far. 
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. 5t 1. Introduction mental parameters (mass, chemical composition of parent star, 

dynamical environment). 

The increasing number of extrasolar planets discovered by ra- 

j. , , .. , ir\A ID »i » 1 onm More than half of the Stars in the solar neighborhood are 

dial velocity surveys (over 100, see e.g. Butler et al. .2002 ° 

Butler et al.|inn3 Udry et al^M reveals a surprising variety ^^"^'"3' ^^stems. Recent discoveries of planets in binary sys- 

of planetary systems: giant planets very close to their parent ^^"^^ 'liff^i"^"^ "^^^^ 1"^^^°^ ^"^^ separations indicate that 

stars, planets in very eccentric orbits, multiple planets systems, Pl^"^^^ "^ig^t form in binaries. The case of 16 Cyg B is re- 

, . , .. 1 c »u J J- t J markable, showing radial velocity variations due to a massive 

in some cases showing orbital resonances. Further dedicated ° ^ 

re . ■ J f ff planet orbiting the star at 1.7 AU (Cochran et al. A991i . The 

efforts are required for a detailed understanding of frequency f ^ ^ i ' 

■ f • . 1 . f ^- f f J very high eccentricity (e=0.68) might be induced by dynam- 

and properties or giant planets as a function of some funda- ^ ^ ^ / & ^ 

ical perturbations of 16 Cyg A (Holman et al. 119971 , Mazeh 

et al. 1997) . Even more compelling is the case of y Cep, that 

Send ojfprint requests to: S. Desidera harbours a planet in spite of the presence of a low-mass stel- 

* Based on observations made with the Italian Telescopio Nazionale lar companion in a 70-year orbit (Hatzes et al. 2003). Zucker & 

Galileo (TNG) operated on the island of La Palma by the Centro Mazeh (2002.) noted that planets found in binaries have a differ- 

Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the gnt mass-period correlation than planets orbiting single stars, 

Spanish Observatorio del Roque de los Muchachos of the Instituto ^-^^ ^„ overabundance of short-period high-mass planets, 
de Astrofisica de Canarias, and observations collected at the European 

Southern Observatory, Chile, using FEROS spectrograph at 1.5 m tele- Besides the study of dynamical effects (minimum binary 

scope (proposal ID: 69.D-0338) separation for planets formation, perturbations of planetary or- 
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bits), binary systems could be used to explore the origin of the 
connection between high metallicity and occurrence of plan- 
etary companions (Gonzalez et al. 120011 Santos et al. 2001 1. 
Two possible explanations were suggested: either a high metal 
content makes it easier to form planets (and in this case no 
metallicity difference is expected between members of a bi- 
nary system with and without planetary companions) or the 
high metallicity is the result of planets or planetesimal inges- 
tion (Gonzalez II 997> . It is also possible that both explanations 
are valid. 

Infall of some planetary material on the star surface is 
expected during the formation and subsequent evolution of a 
planetary system. This may happen during the early phases of 
the planet formation, or in a more mature stage of the planetary 
system evolution, with the infall of remnant of the protoplane- 
tary disk or with the ingestion of already formed planets, scat- 
tered toward the star by dynamical interactions with other plan- 
ets in the system, distant stellar companions or passing stars. 

If the material engulfed by the star is depleted in hydrogen 
and helium, the chemical composition of the star can be modi- 
fied, leading to metal rich convective envelopes. This alteration 
is small or negligible for stars in the pre main sequence phase, 
that have thick convective zones. Stellar pollution is likely effi- 
cient mostly for main sequence stars with masses 1.0-1.3 Mq, 
that have a convective envelope thin enough to allow a signif- 
icant metallicity enhancement for reasonable quantities of ac- 
creted material. 

The search for stellar pollution signatures was prompted by 
the evidence that the host stars of the extrasolar planets found 
by radial velocity surveys are typically metal-rich. The stars 
with planets were searched for trends in metallicity as a func- 
tion of stellar mass or size of convective envelope (Santos et 
al. 1200 II . Pinsonneault et al. 20011, for anomalous abundance 
ratios, signatures of accretion of chemically fractionated ma- 
terial (Smith et al. I2001> . and for Uthium content anomalies 
(Rvan l2000l Israehan et al.'200r Reddy et al.''200Tl. Lithium, 
and particularly the ^Li isotope is in fact a good diagnostic 
for accretion since it is destroyed during the stellar lifetime. 
Murray et al. (2001) studied instead field main sequence stars 
with spectroscopically determined iron abundance. They found 
evidences for typical accretion of 0.5 of iron during the 
main sequence lifetime. 

However, these results are hardly conclusive because of the 
lack of a suitable reference representing the original, unpol- 
luted abundances. 

Wide binaries represent instead easier targets to search for 
chemical anomalies, because of the comparison immediately 
provided by the two companions. Furthermore, their physical 
association (i.e. same distance) allows the use of differential 
analysis between the components. In fact, some of the most 
compelling results on the chemical anomalies of the stars with 
planets were achieved by studying the stars with planets mem- 
ber of binary systems e.g. 16 Cyg (King et al. 119971 Laws & 
Gonzalez'MJB and HD 178911 (Zuckeret al'MHI. 

Differential techniques for the determination of stellar 
abundances are very sensitive for pair of stars with similar tem- 
peratures and gravities, i.e. main sequence stars with similar 
masses. This argument was at the basis of our project for the 



study of a sample of wide binaries, aimed to search for planets 
(using radial velocities) and the chemical composition signa- 
tures due to the engulfment of planetary material. This survey 
is on-going at TNG (the Italian Telescopio Nazionale Galileo, 
INAF, on La Palma, Canary Islands, Spain), using the high res- 
olution spectrograph SARG (Gratton et al. 2001ai. 

The first result of our survey was the discovery of a pair 
(HD 2 1 9542) whose components differ in iron content by about 
0.1 dex (Gratton et al. 2001b hereinafter Paper I), with hints 
that the overabundance of the primary is confined to rocky el- 
ements, while volatiles have comparable abundances. The pri- 
mary is also lithium rich. The most likely explanation of this 
abundance pattern is the infall of planetary material on the pri- 
mary, either portions or remnants of a circumstellar disk or al- 
ready formed planets. Reddy et al. ( I2002I I searched for *Li in 
HD 219542 A. This isotope represents a very powerful trace of 
recent accretion. They find ^Li/^Li=0.03+0.03, clearly not con- 
clusive. It must be noted that such an amount of ''Li is nearly 
the one expected on the basis of the observed iron difference. 

In this paper, we present the results of the 2 years radial 
velocity monitoring of the components of HD 219542 aimed 
at finding planets around them. In this case, we are searching 
for planets around stars showing evidence of accretion of plan- 
etary material. This effort is complementary to the searches for 
abundance anomalies of stars around which planets have al- 
ready been detected, and may represent an important piece in 
the understanding of the origin of the observed chemical com- 
position differences. 

The paper is organized as follows: in Sect. 2 and 3 the avail- 
able information concerning the stellar and orbital characteris- 
tics of the HD 219542 binary system are described; in Sect. 4. 
our high precision radial velocity measurements are presented; 
in Sect. 5 the significance of the variability of HD 219542 B 
is discussed; in Sect. 6 we discuss the possibility that the ob- 
served variations are due to activity-related phenomena; in 
Sect. 7 upper limits to the planets still allowed to exist around 
HD 219542 A are derived; in Sect. 8 the implications of our 
result are discussed. 

2. Stellar properties 

HD 219542 (HIP 114914, BD -02 5917, ADS 16642) is a 
wide pair of main sequence solar type stars. The physical asso- 
ciations of the components was confirmed by the Hipparcos 
astrometry. The main properties of the components of HD 
219542 are summarized in Tabled 

2.1. Stellar masses and ages 

Isochrone fitting can only poorly constrain the age of unevolved 
main sequence stars. The position in the color-magnitude dia- 
gram is best fitted by an age in the range 1-2 Gyr The other 
indirect age indicators Uke chromosperic activity. X-ray emis- 
sion, lithium, and rotational velocity are discussed in Sect. 12.21 
12.312.41 

The stellar masses derived using the Z=0.03 1 and 2 Gyr 
isochrones of Girardi et al. ("MKri are 1.15-1.12 and 1.07-1.05 
for A and B component respectively, with an uncertainty of 



S. Desidera et al.: A search for planets in the metal-enriched binary HD 219542 



3 



Table 1. Stellar properties of the components of HD 219542. 2.2. Rotational profiles 



Parameter 



HD 219542 A HD 219542 B Ref. 



a (2000) 


23 16 35.24 


23 16 35.42 


1 


6 (2000) 


-01 35 11.5 


-01 35 07.0 


1 




163.28 


+ 2.12 


163.28 ±2.12 


1 


fis (mas/yr) 


20.97 H 


= 1.52 


20.97 ± 1.52 


1 


RV (km/s) 


-11.1 : 


t 0.2 


-10.0 ± 0.2 


2 


n (mas) 




18.30 ± 1.97 


1 


rf(pc) 




54.6 ± 5.9 


1 


U (km/s) 




-40.5 d 


:4.8 


2 


V (km/s) 




-16.4 d 


:3.6 


2 


W(km/s) 




-4.0 + 


2.5 


2 


Vt 


8.172 ± 


0.019 


8.591 ± 0.028 


3 


(B - V)t 


0.705 ± 


0.029 


0.791 ± 0.044 


3 


V 


8.198 ± 


0.021 


8.589 ± 0.021 


4 


V-I 


0.641 ± 


0.035 


0.705 ± 0.037 


4 


Hp 


8.312 ± 


0.006 


8.689 ± 0.008 


1 


Hp scatter 




0.01) 


a.h 


1 


J 


6.992 ± 


0.023 


7.216 ± 0.033 


5 


H 


6.726 ± 


0.039 


6.896 ± 0.045 


5 


K 


6.627 ± 


0.033 


6.830 ± 0.045 


5 


b-y 




0.405" ± 


0.006 


6 


mi 




0.214" ± 


0.009 


6 


Cl 




0.367" ± 


0.004 


6 



ST 


G2V 


G7V 


7 


Mv 


4.48 ± 0.23 


4.90 + 0.23 


8 


Teff (K) 


5989 ± 50 


5713 + 50 


8 


^T,ff(A - B) (K) 


276: 


t 14 


8 




4.37 ± 0.10 


4.38 + 0.10 


8 




-5.36 + 0.13 


-5.11 +0.13 


2 


V sin I (km/s) 


1.9 + 0.1 


2.4 ± 0.1 


2 


Vmac (km/s) 


4.4 ± 0.3 


4.8 ± 0.3 


2 


Lx (erg/s) 


2 1+'-'* 


1028 . 


2 


[Fe/H] 


+0.29 + 0.10 


+0.20 + 0.10 


8 


A[Fe/H](A - S) 


0.091 +0.009 


8 


logA^L, 


2.35 


<1.0 


8 



Mass(Mo) 

Radius(i?o) 
Age (Gyr) 



1.13 + 0.05 
1.05 



1.06 + 0.05 
0.96 



References: 1 Hipparcos (ESA I1997> : 2 This Paper: 3 Tvcho (ESA 
[TWTi : 4 Cuypers & Seggewiss (T999J; 5 2MASS; 6 Olsen {1958); 7 
Corbally & Garrison IT980t : 8 Paper I 



" A+B 

^ See Sect.lOlfor details. 



0.05 Mq due to the errors in the parallax. In any case, the 
derivation of ages and masses using standard stellar models can 
be misleading for polluted stars (as likely HD 219542 A), that 
should have different composition in the radiative and convec- 
tive zone (Sasselov i2002i Sandquist et al. i2002ii . 



Our template spectra have resolution {R ~ 150000) and S/N 
(~ 200) high enough that the shape of the Fourier Transform 
of portions of spectra including suitable lines can be used to 
derive rotational and macroturbulent velocities. 

Stellar rotational velocities have then been derived for both 
components of the system (A and B) using the first zeros 
of the Fourier Transform, as described by Gray 019921 To 
this purpose, we used the lines at 6252.565, 6421.360 and 
6439.083 A measured on the radial velocity template spectra. 
Assuming a limb-darkening coefficient of 0.6 (appropriate for 
the stars under consideration), we obtain estimates of 1 .9 + 0.1 
and 2.4 + 0.1 km/s for the values of v sin / for components A 
and B respectively. The value for component A agrees quite 
well with the value of 1.55-2.24 km/s obtained by Reddy et 



al. (I2002> . 

The macroturbolence velocities result of 4.4 + 0.3 and 4.8 + 
0.3 km/s for HD 219542 A and B respectively, normal for this 
type of stars. 

It should be noted that there is some degree of degeneracy 
between the values for rotational and macroturbulent veloci- 
ties, that is fully removed only by using very high S/N spectra. 
However, the macroturbolence we obtain is consistent with the 
spectral types of the stars, supporting the adopted values for 
rotational velocities. 

2.3. Activity indicators 

The level of chromospheric activity can give some clues about 
the age of the system and it is relevant also for the interpreta- 
tion of the radial velocity behaviour of the stars (see Sect. |6j. 
We acquired spectra of both components of HD 219542 (as 
well as of other ~ 40 stars) using the FEROS spectrograph at 
the 1.5 m telescope in La Silla in August 2002. Details of ob- 
servations will be reported elsewhere. The spectra cover the 
range 3900-9200 A, including the region of chromospheric di- 
agnostics H and K Call lines. Fig. [2 shows the region of the 
K line for HD 219542 A and B, compared to the solar spec- 
trum and to the active star 36 Oph A. It appears that both stars 
have low chromospheric activity, with some emission possibly 
present in component B and no evidence for emission in com- 
ponent A. In order to quantify the activity level, we measured 
on our spectra the fluxes within a 1 A band centered on the Ca 
II K line, and normalized them to the values of the fluxes in 
two narrow pseudocontinuum bands at 3909 and 3954 A. This 
procedure was repeated for all stars observed with FEROS, as 
well as for solar type stars observed during the FEROS com- 
missioning (Kaufer et al. .19998 . whose spectra are available 
on the FEROS Consortium archives'. The resulting flux ratios 
I(K) were converted into R^^j, values by using 41 observations 
of 1 1 stars with values of 7?^^ determined in the literature (from 
Baliunas et al. [19951 Henry et al.ri996 and Henry et al. l2000t . 
The transformation included a term dependent on (B-V); this 
term accounts for the systematic dependences of the fluxes in 
the narrow reference bands on temperature and metal abun- 



http;//feros. lsw.uni-heidelberg.de/docu/ferosDB/search.html 
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dance. The finally-adopted correlation is able to predict R^j^ 
values for the calibrators with an accuracy of 0.10 dex (see 
Fig. 2}; the residual scatter is very similar to the r.m.s. scatter 
of R^uj^ values for the Sun taken at different epochs (Baliunas 
et al. [T^95t . and may be attributed to the variability of the Ca II 
emission^. The values of r!^j^ for HD 219542 obtained by this 
procedure are -5.36 + 0.13 and -5.11 + 0.13 for components 
A and B respectively. This low level of chromospheric activity 
is consistent with the results quoted by Glebocki et al. ( 1980 1. 
The small rotational velocities we measured are also in agree- 
ment with a low activity level for both the stars. 

The ROSAT satellite detected a source (IRXS J231636.1- 
013459) that has an offset of 13 and 18 arcsec with re- 
spect HD 219542 A and HD 219542 B respectively (Voges et 
al.|2000l. Since the position error of the X-ray source is 19 arc- 
sec, IRXS J231636. 1-013459 is very likely the counterpart of 
one or both the components of HD 219542 (at the high galactic 
latitude of this star the chance of a background object of similar 
flux is very small, according to Hasinger et al. il993i) . The X-ray 
luminosity derived using the calibration of Hiinsch et al. (I1999> 
is Lx = 2.1;^[ 2 10^*^ erg/s (A+B components). Assuming simi- 
lar contribution by both components, the X-ray luminosity re- 
sults about one order of magnitude lower than Hyades stars and 
typical for stars with an age of about 1-2 Gyr (Randich 2000 1. 
Such X-ray luminosity is also in reasonable agreement with 
the observed rotational velocities according to the relation by 
Bouvier ( flWOt . 

2.4. The galactic orbit and the possible association 
with the Hyades 

Eggen (119601 ) in his search for moving groups in the so- 
lar neighboorood picked up HD 219542 as a member of the 
Hyades moving group. While the existence itself of moving 
groups is debated in the literature (see Taylor 2000 and refer- 
ences therein), we tried to confirm or reject this possible asso- 
ciation using updated velocities and distance for HD 219542 
(see Table^. The space velocities of the binary are (U, V, W) = 
(-40.5, -16.5, -4.0). This is compatible with the space veloci- 
ties of the Hyades given by Ferryman et al. ( 1998 1. The galac- 
tic orbit of HD 219542, computed as in Barbieri & Gratton 
J2002> . is represented in Fig|3] The orbit (Rmin = 7.70 kpc; 
Rnwx = 9.41 kpc; Z,m.i = 0.06 kpc; e = 0.10) is normal for a 
moderately young star. 

The distance from the Sun of HD 219542 (54.6 pc) is com- 
parable to that of the Hyades. The physical distance of HD 
219542 from the center of the Hyades is about 40 pc. 

Other physical indicators can be used to shed light on the 
possible link between the HD 219542 system and the Hyades: 

- The metallicity of HD 219542 B (assumed to be unpol- 
luted) is [Fe/H] =-1-0.20 + 0.10, compatible within errors 

^ Actually, we expect that our values of i?^^ contain systematic er- 
rors due to stellar rotation. Lines of some of our calibrators are indeed 
significantly broadened by rotation; however the values of R'(HK) we 
derived for these stars do not differ systematically by more than 0.05 
dex from the average calibration of Figure|2| 



1.5 



0.5 









1 1 

■ 

h ! 








r 

HD219542B /U 






[ 1 



J930 39.32 3934 

Wovelength (A) 



Fig. 1. The region near the core of the Call K line for the Sun 
(logR'jj^ - -4.94), the components of HD 219542 and the ac- 
tive star 36 Oph A (log 7?^^ = -4.58, Henry et al. .l996j . The 
spectra for the other stars were arbitrarily shifted in intensity 
with respect to the Solar one. The horizontal lines represent the 
zero level for each spectrum. The vertical lines show the spec- 
tral region used for chromospheric activity measurement (the 
continuum regions are outside the plotted range). 

with the cluster metallicity ([Fe/H]=-i-0.14, Ferryman et 
al. [T^ 3 

' Note that the [Fe/H] values given in Paper I may be slightly over- 
estimated, due to the choice of a very low value for the microturbulent 
velocity. 




log R'(HK) (lit) 



Fig. 2. The calibration of R'^^ index based on 1 1 stars with lit- 
erature determinations, observed with the FEROS spectrograph 
during our observing run (August 2002) and instrument com- 
missioning (Kaufer et al. .l999j . 
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r (kpc) 



Fig. 3. Galactic Orbit of HD 219542 in the meridional plane. 



- There are no peculiar abundance patterns either in the 
Hyades stars or in HD 219542 B 

- The hthium content of HD 219542 and especially of 
HD 219542 B is lower than that of Hyades stars of similar 
temperature (Fig. |4j and rather similar to that of stars in 
older open clusters like NGC 752 and M67 (Paper I). This 
seems not compatible with the Hyades association, since 
Hyades show a very low scatter in lithium abundance at a 
given temperature. 

- Isochrone fitting does not have the accuracy for distentan- 
gling ages differences of ~ 0.5 Gyr for unevolved objects 
like the two components of HD 219542 

- The rotational velocities of the components of HD 219542 
are at the low end of those of Hyades stars of similar spec- 
tral type. Note that a pole-on orientation for both compo- 
nents is unlikely, since for binaries with separation larger 
than 40 AU the rotation axes do not appear to be aligned 
(Hale[T994). 

- The chromospheric activity as measured on our FEROS 
spectra is much lower than that of Hyades stars of similar 
colour 

- The X-ray luminosity derived in Sect. l2.3l is about one or- 
der of magnitude lower than that of Hyades stars of similar 
temperature 

- None of the 24 Hyades G dwarfs studied by Paulson et 



3.5 
3.0 



1.5 
1.0 h 
0.5 
0.0 



8 8 



V- 

^ I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

6400 6200 6000 5800 5600 5400 5200 5000 

Tef, (K) 



Fig. 4. Lithium determinations for Hyades stars from Thorburn 
et al. f 1993'! (empty circles for abundance measurements and 
empty triangles for upper limits) and for the components of 
HD 219542 (HD 219542 A: filled circle; HD 219542 B: filled 
triangle). Only an upper limit was derived for HD 219542 B. 
The line at the lower left corner represents variations of derived 
lithium content as a function of the adopted effective tempera- 
ture. Even assuming large differences in the temperature scales 
used by Thorburn et al. (1993 1 and Paper I of ~ 150 - 200 K, 
it is not possible to obtain a lithium content for HD 219542 
similar to that of the stars of Hyades of similar temperature. 

Table 2. Relative position of the components of HD 219542. 



Epoch 


P 


Error 


e 


Error 


^^2000.0 


Ref. 


1830.51 


4.56 




26.7 




26.5 


1 


1890.64 


4.75 




30.6 




30.5 


1 


1909.05 


4.96 




29.9 




29.8 


1 


1917.32 


4.96 




30.3 




30.2 


1 


1924.40 


4.97 




29.2 




29.1 


1 


1953.545 


5.103 


0.013 


30.94 


0.16 


30.89 


2 


1960.699 


5.11 


0.02 


30.8 


0.2 


30.8 


3 


1976.784 


5.02 




31.3 




31.3 


4 


1991.25 


5.279 


0.004 


31.53 




31.52 


5 


1991.797 


5.280 


0.002 


31.47 


0.02 


31.46 


6 


1994.747 


5.29 


0.20 


32.05 


0.78 


32.05 


7 



al. J2002> has a radial velocity jitter as low as HD 219542 A 
(see Sect.lOli. 



References: 1 Aitken a932j ; 2 Van Albada - Van Dien (1978); 3 Van 
Albada - Van Dien (M^}; 4 Holden {WT\ : 5 Hipparcos (ESA fTWl 
6 Cuypers & Seggewiss <1999> : 7 Popovic & Pavlovic <1997t . 



In summary, we conclude that while the kinematic data 
suggests the HD 219542 may be part of the Hyades moving 
group, the other stellar properties (lithium, rotation, chromo- 
spheric activity. X-ray emission, radial velocity jitter) seems to 
indicate a much older age (at least 1-2 Gyr) for the system. 
In particular, isochrone fitting and X-ray emission point to an 
age of ~ 1 - 2 Gyr, while the very low chromospheric activity 
suggests an even older age. 

The measurement of lithium of HD 219542 A of Paper I was con- 
firmed by the independent analysis by Reddy et al. c2002< . 



3. Orbital properties of tlie binary system 

3.1. Astrometric Data 

The projected separation of 5.279 arcsec at the distance mea- 
sured by Hipparcos (54.6 pc) corresponds to 288 AU. Available 
data can be used to costrain the orbit of the binary. The first as- 
trometric measurement of this pair dates back to 1830 (Aitken 
I1932> . The literature astrometric measurements are listed in 
Table |2l Position angles were precessed to 2000.0 using Eq. 
5ofHeintz(.1978.). 
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Fig. 5. Astrometric measurements of HD 219542. Overplotted 
is the best linear fit to the data. 

Clear trends both in separation and in position angle can be 
noticed (Fig. |5|i. No changes in separation and position angle 
were instead detected by Hipparcos because of its limited time 
baseline. A linear fitting to the available data^ gives dpjdt - 
0.0048 + 0.0003 arcsec/yr and Mjdt = 0.018 + 0.003 deg/yr. 

3.2. Radial Velocity Difference 

A further clue to the orbital properties of the system 
comes from the radial velocity differences between the 
components, that results -1.15 + 0.02 km/s (see Sect. 4). 
Therefore HD 219542 B is moving toward us with respect to 
HD 2 19542 A. 

3.3. Clues on the orbital parameters 

The availability of p, 6, their time derivatives, the velocity dif- 
ference and the parallax means that we know the relative posi- 
tion of the stars in the plane of the sky (x, y) and all the velocity 
components Vx, Vy, v, (Table|3}- The position z along the line of 
sight is instead unknow, so that it is not possible to derive the 
orbital parameters of the visual binary. 

However, it is possible to constrain the value of z assum- 
ing that the two components are bound and derive the orbital 
parameters as a function of z- To this aim, we follow the proce- 
dure by Hauser & Marcy (I1999> . The allowed separation along 
the line of sight results -1260 < z < 1260 AU. The orbital 
parameters as a function of z are shown in Fig.|6l 

It appears that low eccentricity orbits with periastron of a 
few hundreds AU, as well as very eccentric orbits with perias- 
tron as small as 10 AU, are compatible with the available data. 
Radial velocity variations due to the binary orbit are below 1 
m/s/yr over the whole range of z, undetectable from our data. 

The critical semiaxis for dynamical stability of planets for 
the family of orbits described above can be calculated follow- 

^ The measurements were weighted according to their errors, when 
available, and assigning errors of 0.2 arcsec and 2 deg for p and 6 
measurements respectively for those stars without errors bars. 



Table 3. Relative position and velocities of the components of 
HD 219542 in the cartesian plane. 



Parameter 


Value 


Error 


x(AU) 


246 


26 


y(AU) 


151 


16 


Zmln (AU) 


-1260 


280 


Z^x (AU) 


1260 


280 


Vx (m/s) 


840 


140 


Vy (m/s) 


1020 


110 


Vz (ni/s) 


-1150 


25 




z (AU) z (AU) 

Fig. 6. Possible values of semimajor axis, period, eccentricity 
and periastron distance of the orbit of HD 219542 as a function 
of the (unknown) separation along the line of sight. 

ing Holman & Weigert C1999) . As expected, for the very ec- 
centric orbits occurring for -500 < z < AU the zone allowed 
for stable planetary orbits is very small, while for the low ec- 
centricy orbits occurring if z ~ 500 AU dynamical stability up 
to ~ 100 AU is ensured (Fig.|3. 

4. Radial Velocities 

4.1. Observations and data analysis 

The components of HD 219542 were observed as part of our 
ongoing planet search program using S ARG at TNG. The goals 
of the project are the study of the frequency of planets around 
stars in wide binaries with atypical separation 200-300 AU, 
looking for the dynamical effects due to the stellar compan- 
ion, and the search for chemical differences between the com- 
ponents (see Paper I). A recent report is presented by Gratton 
et al. (2002). The survey as well as the observational and data 
analysis procedure will be described in detail in a forthcoming 
paper (Desidera et al., in preparation). Here we briefly summa- 
rize the most relevant points. 

The observations used for the radial velocity determi- 
nations were acquired with the SARG spectrograph using 
the Yellow Grism, that covers the spectral range 4600- 
7900 A without gaps, and the 0.25 arcsec slit. The resulting res- 
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Fig. 7. Critical semiaxis for dynamical stability of planets or- 
biting HD 219542 A (upper panel) and HD 219542 B (lower 
panel) as a function of the current separation along the hne of 
sight. 

olution is R=150000 (2 pixels sampling). The iodine cell was 
inserted into the optical path, superimposing a dense forest of 
absorption lines used as reference spectrum for the radial veloc- 
ity determination. The spectral ranges 4600-5000 A and 6200- 
7900 A are free from the contamination of iodine lines and can 
be used to monitor the changes of the line profile and equiva- 
lent width (see Sect l6.1l - |6.3> . Only the blue part of the spectra 
is used in this paper The insertion of a suitable filter avoids 
any grism second order contamination in the red extreme of 
the spectra. Exposure times were fixed at 15 minutes, to re- 
duce the errors in barycentric correction caused by the lack of 
knowledge of the exact flux mid time of the exposure. During 
the observations, the slit was usually oriented perpendicularly 
to the separation of the components to minimize the contami- 
nation of the spectra by the companion. The rather large sepa- 
ration of the components of HD 219542 (5.3 arcsec) makes the 
contamination issue not important for the case presented here. 

The data reduction was performed using IRAF^. 

The radial velocity analysis was performed using the 
AUSTRAL code developed by Endl et al. (.2000.) . Each 
starH-iodine spectrum was modeled using a pure iodine tem- 
plate of the SARG cell obtained with the Kitt Peak FTS 
(Desidera 2000 1 and a high quality pure stellar template (i.e. 
without the iodine cell in the path) obtained with SARG, prop- 
erly deconvolved for the instrument profile. In the modeling of 
the StarH-iodine spectrum, the best Doppler shift of the stellar 
spectrum relative to its template is found and the instrument 
profile was reconstructed. Since the instrument profile changes 
along a spectral order, the spectrum was divided in chunks ap- 
proximately 100 pixels long (~2 A), each one modeled inde- 
pendently. About 400 chunks were used in the final analysis. 
Internal errors of velocities are calculated as the errors of the 
mean of the chunks. We find that averaging the velocities ob- 

* IRAF is distributed by the National Optical Observatory, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Fundation 
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Fig. 8. Differential radial velocities of HD 219542 A. The dis- 
persion of the velocities is 5.8 m/s, compatible with measure- 
ment errors. 

tained with several different chunk lenghts (in the range 90-120 
pixels for solar type stars) improves the quality of the results. 
Changing the chunk length means to slightly reduce the impact 
of noise due to the division of the spectrum in discrete pieces. 

The parameters for the instrument profile modeling were 
instead kept fixed on values fine-tuned using fast rotating fea- 
tureless stars, the constant radial velocity star r Ceti and the 
planet-bearing stars 5 1 Peg and p CrB. A long term precision of 
about 3 m/s was achieved for these stars. For the much fainter 
(V=8-9) stars of the binary sample the errors are typically about 
5 m/s, due to the lower S/N of the spectra. For the components 
of HD 2 1 9542, the typical S/N we achieved is 80- 1 00. The stel- 
lar templates (used also in Paper I for the abundance analysis) 
have instead S/N above 200. 

4.2. Differential Radial Velocities 

More than 20 radial velocities measurements of the com- 
ponents of HD 219542 were obtained, spanning the period 
October 2000 - November 2002 (TablesgE}. On a few nights, 
we obtained two consecutive spectra for the targets. In these 
cases, we computed a weighted average of the velocities. 

The time series of velocities of HD 219542 A is shown in 
Fig. |8l The dispersion of velocities is 5.8 m/s, only slightly 
larger than the mean internal errors (5.6 m/s) determined by 
the chunk-to-chunk internal scatter in the velocities^. Statistical 
analysis (see Sect. |5} does not reveal significant periodicities. 
Long term trends are also not apparent from the data. Therefore 
we can rule out both short-period Jupiter-mass planets as well 
as massive companions in external orbits. The upper limits on 
the giant planets around HD 219542 A still compatible with the 
data are discussed in Sect.0 

' Comparison between pairs of observations obtained consecutively 
in both components suggest that the errors given by this procedure 
may be overestimated by about 25%, possibly because individual 
chunks consistently deviates by similar amounts in spectra acquired 
on the same night. 
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Table 4. Differential radial velocities of HD 219542 A 



JU -z4jUUUU 


T>\7 
KV 


error 


1825.52 


-7.2 


5.8 


1826.59 


-6.0 


8.6 


2070.70 


-1.9 


3.8 


2071.70 


1.9 


3.9 


2072.68 


1.0 


4.5 


2113.69 


-2.4 


4.1 


2115.67 


-3.8 


4.5 


2115.71 


4.0 


5.1 


2116.56 


5.8 


8.2 


2116.70 


-1.4 


5.6 


2117.68 


5.2 


3.9 


2120.71 


1.1 


5.0 


2145.65 


-1.4 


4.0 


2190.53 


4.0 


4.5 


2216.46 


-5.5 


7.7 


2245.42 


-3.1 


4.4 


2424.70 


3.7 


4.6 


2445.70 


11.6 


4.5 


2472.70 


7.9 


6.0 


2504.67 


-7.2 


11.0 


2538.50 


-7.5 


5.9 


2570.40 


-2.1 


4.9 


2585.46 


12.9 


7.3 


2597.36 


-2.3 


5.2 


2605.35 


-7.3 


5.1 



The velocities of HD 219542 B show instead a dispersion 
(9.2 m/s) larger than the measurement errors (5.1 m/s). Since 
the two components were usually observed consecutively, we 
can exclude the occurrence of systematic errors due to uncor- 
rected instrumental shifts or barycentric correction as possible 
causes of the observed variations. 

We performed a period search within the RV-results in the 
range of 2 to 1000 days, using the well-known Lomb-Scargle 
periodogram (Lomb 1976 Scargle '1982 i. We found a strong 
peak at a period of 112 days. Its significance is discussed in 
Sect. |5] The best-fit Keplerian orbital solution^ to the RV-data 
of HD 219542 B yields a period P=l 12. 1 days, semiamplitude 
K=13.0 m/s and a moderate eccentricity e ~ 0.32 (Table|6j. 
The r.m.s. of residuals from the keplerian fit is 5. 1 m/s (5.9 m/s 
taking into account the loss of degrees of freedom introduced 
by the keplerian fitting). Fig. |5] shows the radial velocities of 
HD 219542 B with overplotted the orbital solution; Fig. [TOI 
shows the residuals from the keplerian fit, and Fig.[nthe ve- 
locities of HD 219542 B phased to best fitting period. Fischer 
et al. (1200 U noted that for noisy data (and in the case of un- 
detected further planets in the system) the orbital fitting pro- 
duces spuriuos eccentricity enhancements. For this reason, the 
derived eccentricity should be considered with some caution. 

Among the planets found so far, 3 have comparable veloc- 
ity amplitude (namely HD 16141b, Marcy et al. .2000. 47 UMa 
c, Fischer et al. l2002l 55 Cnc c, Marcy et al. l2002t . Therefore 

* The fitting was performed using the program Orbitsolver (R.E.M. 
Griffin, priv. comm.) that adjusts all the elements of the orbit simulta- 
neously in a least-squares solution (Griffin & Comell il99^ . 
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Fig. 9. Differential radial velocities of HD 219542 B. 
Overplotted is the best keplerian solution. 
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Fig. 10. Residuals of the velocities of HD 219542 B from the 
calculated orbit. 
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Fig. 11. The velocities of HD 219542 B phased to the 112.1 
days period of the best orbital solution. 
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Table 5. Differential radial velocities of HD 219542 B 



JJJ -z4jUUUU 


KV 


error 


1825.51 


-4.4 


7.5 


1826.49 


-11.9 


5.9 


2070.71 


1.5 


4.1 


2071.71 


6.6 


4.0 


2072.69 


8.9 


5.9 


2113.71 


-0.2 


4.6 


2115.69 


2.4 


5.7 


2115.71 


8.6 


5.2 


2116.68 


-6.6 


5.8 


2117.70 


-3.8 


7.8 


2117.71 


2.5 


6.8 


2120.72 


2.4 


5.0 


2145.66 


-2.2 


4.6 


2190.54 


3.1 


3.8 


2216.47 


9.4 


7.2 


2245.43 


-8.4 


4.3 


2424.71 


11.4 


5.2 


2445.71 


13.2 


4.5 


2472.70 


-8.6 


5.2 


2538.51 


23.4 


7.9 


2570.41 


-7.8 


5.2 


2570.42 


-0.7 


4.9 


2585.44 


-16.4 


7.3 


2585.45 


-18.4 


6.7 


2597.33 


-1.9 


4.4 


2597.34 


5.3 


4.8 


2605.36 


-4.8 


4.8 


2605.37 


-2.8 


5.3 



Table 6. Possible Orbital Solution for the velocity of 
HD 219542 B. 



Parameter 


Value 


Period (d) 


112.1 + 1.2 


Semi-amplitude (m/s) 


13.0 + 3.0 


Eccentricity 


0.32 ±0.15 


Longitude of periastron (deg) 


51±33 


Periastron passage (JD-2450000) 


1770.6 + 8.0 


a(AU) 


0.46 


m sin ; 


0.30 



a careful evaluation of the data is required to discuss the sig- 
nificance of the result (Sect.|5} and the possible activity origin 
of the observed variations (Sect.|6j. As we will see, the results 
of these tests are rather inconclusive, so that the existence of a 
planet around HD 219542 B still needs confirmation^. 

4.3. Absolute Radial Velocities 

The iodine cell technique provides velocities relative to an arbi- 
trary zero point, different for each star. Cross-correlation with 

' We cannot then at present confirm at a high level of confidence 
the announcement of a planet that some of us (the Italian part of the 
group) gave as a press release in November 2002. 



a template of a star with known absolute radial velocity allows 
to determine the absolute velocity values for both stars. 

We used as template 5 1 Peg, that has an absolute velocity 
of -33.23 km/s (Nidever et al. 2002 1. The velocities result of 
Va - -11.1 km/s and Vb - -10.0 km/s. Corrections for grav- 
itational redshift and convective blueshift were not included, 
because 5 1 Peg has a colour intermediate between those of HD 
219542 A and B. It was also not necessary to include a cor- 
rection for the orbital motion due the planet around 51 Peg, 
because the template of 5 1 Peg was acquired by chance with 
the star close to its center of mass velocity. The eiTors of our 
velocities are likely within 0.2 km/s. The only literature de- 
termination of radial velocity is -16.5 km/s for HD 219542 A 
(Wilson ri953> . The discrepancy is likely explained by the low 
accuracy of the old data. 

For the analysis of the possible binary orbit (Sect. |3Jl it is 
important to have a precise measurement of the difference of 
radial velocities between the components. Such a difference 
was measured using the same stellar template for both stars. 
Considering the templates of HD 219542 A and HD 219542 B, 
the radial velocity difference results Va-Vb = -1200+10 m/s. 
Such quantity must be connected for the difference of the grav- 
itational redshift and convective blueshift between the stars to 
give the real radial velocity difference. We estimate this effect 
to be about 50 m/s, using the relation by Nidever et al. ( I2002> . 
Therefore our best estimate for the radial velocity difference 
between the components is Va - Vb - -1150 + 25 m/s (in- 
cluding further ~ 20 m/s of uncertainty in the gravitational and 
convective corrections). 

5. Test of significance of tlie proposed planet orbit 

We considered two different tests to evaluate the significance of 
the solution. The first one answers to the question: "What is the 
probabiUty that we would extract a solution corresponding to a 
Keplerian motion induced by a planet better than the solution 
found from the real data from a similar analysis of a random 
sequence of radial velocities sampled at the same epochs and 
with the same noise as the data ?" 

The following procedure was used to answer this question. 
We employed an orbit finder program that scans through the pa- 
rameter range for period, amplitude, phase angle, eccentricity 
and periastron angles and determines the quality of the fit for 
each test-orbit. We considered two different best fit estimators, 
minimizing either the r.m.s. scatter of data or the;\f^ of the dis- 
tribution, using the internal eiTors as estimate of the real errors. 
While might be expected to be a better fit estimator, we note 
that on one side the distribution of individual internal errors 
is quite uniform, and on the other internal errors may neglect 
some source of scatter present in the data (e.g. activity-induced 
jitter, barycentric correction, etc.). Furthermore, we have some 
indications that internal errors are overestimated by ~ 25% 
(Sect. l4~2t . We define as acceptable those solutions whose ec- 
centricity is not larger than 0.8, and there is no phase interval 
longer than 0.25 without any data. 

We verified that our orbit finder program delivers similar 
solutions to the OrbitSolver program using both estimators. 
The solutions found by our program have in both cases period 
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Fig. 12. Lomb-Scargle periodogram results for the RV data of 
component A (upper panel) and component B (lower panel). 
The false-alarm-probability (FAP) levels of 1% are indicated 
as horizontal dashed lines (see text for details). While there is 
no signal present for HD 219542 A, we find a strong peak at 
P= 1 1 1 days for HD 219542 B. However, even this signal has 
a FAP of 3.7% and is thus not highly significant. 

of 111.3 d, amplitude of 12.7 m/s, and eccentricity 0.4. The 
r.m.s. of residuals is of 5.8 m/s and the;^'^ is 104.7 m/s. In both 
cases they are very close to the optimized value. 

We employ a bootstrap randomization scheme in order 
to estimate the significance of the 112 days period. In this 
bootstrap method the radial velocities and errors of individual 
nights are randomly re-distributed along the times of observa- 
tions. 

We ran the program on the generated data sets and we found 
that only in 2 cases out of 1000, random extractions yield so- 
lutions with an rm.s. smaller than or equal to that given by the 
real observations. When using the solutions better than the 
observed case are 1.2%. 

Second, we tested the probability to have by chance a power 
in the periodogram similar to that of the observed data, follow- 
ing the bootstrap approach decribed above and computing the 
power spectrum for each new dataset is generated (Kiirster et 
al. ll997l Murdoch et al. '19931. For each of the two components 
we performed 10000 bootstrap runs. The results are displayed 
in Fig. El We found a false-alarm-probability (FAP) of 3.7%. 
We note that this second approach is more general, since it as- 
sumes that the signal is given by sine-functions, independently 
of its physical origin. From this more prudent estimate, we con- 
clude that there is a 3% probability that our solution may be due 
to chance. 

The same tests performed for component A give a rather 
high probability (more than 50%) that a better result can be 
obtained from a random set of radial velocities. 

Finally, it should also be remembered that a basic assump- 
tion throughout this procedure is that there is no correlation be- 
tween noise on consecutive (or close) nights. However, some 
correlation is expected if noise is dominated by activity. Hence 
this test cannot exclude that the observed velocity variations 
are due to activity rather than to keplerian motion. 



6. Activity or Planet? 

The low amplitude of the velocity variations of HD 219542 B 
may be in principle due to orbital motion (i.e. a planet) or to 
activity related phenomena. We note that a period of 1 12 days 
is too long for a rotational modulation (the upper limit to the 
rotational period from the observed vsin/ value is ~ 21 days) 
and too short for the long term activity cycle typical for solar 
type stars. However, it is still possible that the observed vari- 
ations are due to some transient short period modulations, not 
apparent as significant in the statistical analysis because of the 
lack of coherence over our ~ 800 days of observations In this 
case, however, we would have picked up a phase of low activ- 
ity during our observations in July 2001, when we obtained 
5 observations over 8 nights, showing a dispersion of 4.5 m/s 
(nightly averages). 

We performed a few tests to check this possibility: we stud- 
ied the characteristics of the components of HD 219542 to infer 
the expected activity-induced radial velocity jitter, we checked 
the Hipparcos photometry for variability, and we measured line 
bisector and equivalent width variability on SARG spectra. 

6. 1 . Bisector analysis 

Given the small amplitude of the observed variations in the ra- 
dial velocity of HD 219542 B, it is well possible that such vari- 
ations are not real, but rather are due to the effects of line asym- 
metries. Such variations could be in principle related to activ- 
ity: e.g. in the case of HD 166435, Queloz et al. ( f200H found 
that radial velocity variations with semiamplitude as large as 
A' = 83 m/s can be attributed to correlated variations of the line 
bisectors of ~ 200 m/s (peak-to-peak) with the rotational period 
and attributable to the effects of star spots. Unfortunately, such 
an analysis is much more difficult in the present case, because 
the variations in the bisectors required to explain the radial ve- 
locity variations are very small (< 10 m/s). While we did not 
expect that such an accuracy can be achieved using our data 
for both components of HD 219542, we however performed an 
analysis of the variations of line bisectors, using those orders 
shortward of 5000 A, where there are not lines due to iodine. 

The technique we used is similar to that devised by Queloz 
et al. (2001), that is an analysis of asymmetries in the Cross 
Correlation Function (CCF). The CCF is obtained by cross 
correlating the observed spectra with a suitable mask, that is 
a function that is 1 within very narrow regions (ideally 6 of 
Dirac) centered at the wavelengths of selected spectral lines, 
and elsewhere. This CCF represents an average profile for 
the lines in the lists. To this purpose we selected about 80 clean 
spectral lines in this spectral region, with equivalent widths in 
the range 10-70 mA. The left panel of Fig.^jdisplays a typ- 
ical CCF obtained by this technique. Typical Full Width Half 
Maximum (FWHM) of our CCF's are 0.10 A: this broaden- 
ing is essentially due to macroturbulence, since the instrumen- 
tal profile has a typical FWHM of 0.032 A. In analogy with 
Queloz et al. J200H . we quantified line asymmetry span as the 



The second peak in Fig. ll2l is at 18.5 days (a reasonable rotational 
period for the star). While the season 2002 data can be fitted fairly 
well by a sinusoidal variation with this period, the 2001 data cannot. 
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difference b between the average bisector values in two regions 
of the CCF's (from 0.96 to 0.89, and from 0.86 to 0.72 of the 
Hne profile, see Fig.ll3>. 

Our analysis is certainly not optimal: in fact, asymmetries 
in our CCFs may arise for reasons different from real asymme- 
tries in the spectral lines (and other than photon noise). Possible 
causes include blends with nearby lines, errors in the wave- 
length calibration and inhomogeneities in the slit illumination. 
The first effect has not a large impact here, because it would 
show up as a constant pattern, and we are rather interested in 
bisector variations. Also, we selected lines very carefully, and 
the average line profiles is very close to the continuum at the 
base of the line. We also note that at least in the case of 51 
Peg (for which we performed a similar analysis using our own 
spectra, see below), the run of the bisector with depth is similar 
to those published in the literature (Gray 1 997) in the regions 
of overlap. The runs of the bisectors with depth are very similar 
for 5 1 Peg and the two components of HD 219542, with a trend 
for the span to slightly larger in the warmer star (HD 2 1 9524 A) 
and slightly smaller in the cooler one (HD 219542 B), with 
5 1 Peg (which has intermediate temperature) showing an inter- 
mediate behaviour This agrees with expectations (Grav 1992j . 
The other effects are of major concern, and likely are more im- 
portant in our analysis than in that by Queloz et al. (2001i be- 
cause slit illumination is not scrambled by fibers, and because 
our wavelength calibration is not as accurate. On the other side, 
it should be also noted that our spectra have a resolution that 
is about three times larger and this would result in a gain af 
about one order of magnitude of the sensitivity of the CCF to 
real changes in the bisectors. We also note that all these effects 
do not affect our radial velocities, because they are virtually 
eliminated in our calibration procedure based on the iodine cell 
lines; however they may dominate asymmetries in the cross 
correlation peaks, hiding the effects of activity-related varia- 
tions of the bisectors. 

In order to quantify the internal errors, we analyzed with 
the same procedure a series of 23 spectra of 51 Peg (these 
spectra have S/N ratios about 2.5 times larger than those 
of HD 219542) taken over two years with SARG at TNG. 
Spectrum-to-spectrum variations of the bisector span b are 
~ 10.7 m/s (r.m.s.). The scatter is smaller (about 5.6 m/s rm.s.) 
for spectra taken during the same night (generally over a few 
minutes). This last value seems then a lower limit for the in- 
ternal errors in our bisector values for this star In fact, one 
one hand we expect that real variations of the bisectors should 
show up over longer timescales, and on the other hand we can 
not exclude that other errors arise when comparing data from 
different nights, when e.g. a different wavelength calibration is 
used. Hence, in the case of 5 1 Peg we may set an upper limit to 
real bisector span variations of about 9 m/s from our data. More 
stringent upper limit on the bisector span variations for 5 1 Peg 
was obtained by Hatzes et al. 1999 (~ 7 m/s). However, the 
definition of the bisector span is different in that case, so that 
it is not easy to compare the two results. Transformation of our 
upper limit on bisector to a limit on the jitter in radial veloc- 
ity is not obvious (it depends on the definition of the bisector 
span): however, in the case of HD 166435 (for which Queloz 
et al. 12001 1 used a definition similar to the present one) radial 
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Fig. 13. Left panel: Cross Correlation Function for a spec- 
trum of HD 219542 B; right panel: superposition of the the 
mean line bisectors from the spectra of 51 Peg (dotted line), 
HD 219542 A (dashed Hne), and HD 219542 B (solid line). T 
and B regions are those used in the computation of the bisector 
span. 



velocity variations are roughly 1 . 1 times the variations of the 
bisector spans. If the same rule then applies to our case, the up- 
per limit to radial velocity jitter for 51 Peg is ~ 10 m/s: this is 
about 2.5 times larger than the limit obtained from photometric 
(Henry et al. l2000> and equivalent width variations (see below). 
However, we consider it as still an interesting upper limit. 

In the case of HD 219542, the scatter of bisector spans is 
31.2 m/s and 21.8 m/s for A and B component respectively. 
Bisector span is not significantly correlated with radial veloci- 
ties for both stars (FigO. Lower limits to measurement errors, 
determined as above from the results of consecutive spectra, 
is about 17.3 m/s. Therefore the observed pattern can be ex- 
plained mostly by measurement errors. Anyway, there is still 
room for bisector span variations corresponding to a radial ve- 
locity variations of ~ 15 m/s. This is larger than the observed 
scatter for HD 219542 B, so that we can not exclude that the 
observed radial velocity variations are due to line asymmetries. 

We note that while the bisector analysis is not conclusive 
in the present case, however it would be sufficient to disentan- 
gle the origin of radial velocity variations with amplitudes of a 
few tens of m/s. Finally, we note that a more detailed approach 
(that takes into account the instrumental line profile determined 
using the iodine lines) could be devised; however it requires a 
lengthy procedure and we postpone it to a forthcoming paper. 
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Fig. 14. Bisector span measurements vs radial velocities for 
HD 219542 B (upper panel) and HD 219542 A (lower panel). 
There are not significant correlations for both stars. On the right 
comer the typical short term error on bisector span measure- 
ment is shown. 

6.2. Limits on photometric variations 

The only photometric data set available suitable for a variabil- 
ity search is the Hipparcos photometry. However, it refers to 
the combined light of A and B components. When averaging 
data on a daily basis, the r.m.s. is 0.004 magnitudes, without 
hints of variability above the measurement errors. No signifi- 
cant periodicities appeared applying the same procedures used 
in Sect. |5l The study of the Hyades by Paulson et al. (quoted 



by Henry et al. I2002t gives a relation between the radial ve- 
locity and photometric scatters of ~ 3.3 m/s/mmag. Therefore, 
to justify the observed radial velocity dispersion, a photometric 
scatter 0.003 mag would be required. This is much lower than 
the limit placed by Hipparcos photometry, taking into account 
that HD 219542 B contributes only ~ 40% of the composite 
light. 

6.3. Equivalent Widtli Variations 

We expect that activity may cause (small) variations in the av- 
erage temperature in the line forming region. For this reason, 
we expect (small) variations in the average equivalent widths 
(EWs) of spectral lines. This method is similar to the line depth 
ratio method considered by Gray (1994i. We then selected 
among the line list used in Sect. 16.11 in the spectrum section 
clean of iodine lines the lines with large sensitivity to effec- 
tive temperature changes, determined using appropriate spec- 
tral synthesis. The final selection includes about 40 lines for 
which variations of 1.0 mA in the EW are due to temperature 
changes of less than 20 K. A variations of 20 K in the stellar 
effective temperature causes a variation of about 0.017 mag in 
stellar magnitude. We then measured the EWs of the lines on 
all spectra of both components of HD 219542, as well as for 
the comparison star 5 1 Peg, using an our own purposely writ- 
ten code. On 51 Peg we found that for a single line, the typical 
r.m.s. scatter of the residuals of the EWs with respect of the av- 



erage value is 0.86 mA. The spectra have a typical S/N of about 
170 at this wavelength. Application of the formula by Cayrel 
( 1988) would produce an expected error of 0.22 mA for these 
lines; however this formula neglects the contribution due to the 
continuum and background subtraction, that likely significantly 
increases the errors. Furthermore, internal errors obtained by 
considering the error for individual lines divided by the square 
root of the number of lines, likely underestimates the real error 
bar because they do not take into account correlation existing 
between nearby lines (due to positioning of the continuum level 
and errors in the background subtraction). We then measured 
the temperature variations properly weighting the lines accord- 
ing to their temperature sensitivity. The observed rm.s. scatter 
of results for individual spectra is 2.5 K, with internal errors 
2.1 K. This yields an upper limit to real spectrum-to-spectrum 
of 1 .3 K. This on turn yields an upper limit to real variations 
in the luminosity of 51 Peg of 0.001 1 mag (comparable to the 
results by Henry et al. 120001 . This very strict upper limits con- 
straints the jitter in radial velocity to be less than 3.6 m/s (using 
the quoted relation 3.3 m/s/mmag). 

When the same procedure is applied to the spectra of 
HD 219542 B, we get much less stringent upper limits, due to 
the lower S/N of the spectra. The observed r.m.s. of tempera- 
ture is 8.34 K; the internal error is 5.95 K. This implies an upper 
limit to real variations in the effective temperature of 5.8 K, and 
an upper limit to photometric variations of 0.0050 mag, about 
twice lower than the limit determined from the Hipparcos pho- 
tometry. However, this limit is still large enough to accomo- 
date a jitter in the radial velocity of about 15 m/s, similar to the 
limit obtained from bisectors and well enough to justify obser- 
vations. 

We conclude that also this estimate of the expected jitter is 
not accurate enough for the present purpose. 

6.4. Expected radial velocity jitter 

The low chromospheric activity and the low rotational veloci- 
ties (Rhk = -5.36 + 0.13, vsin/ = 1.9 km/s for HD 219542 
A; Rhk = -5.11 +0.13, vsin / = 2.4 km/s for HD 219542 B) 
indicate a low activity level. The jitter expected on the basis of 
the relation cry vs vsin / by Saar et al. ( 1998) is about 6 m/s 
and 8 m/s for HD 219542 A and B respectively. The logR^^^ 
calibration gives instead ~ 2 and 3 m/s respectively. The cal- 
ibration by Marcy (Marcy 2002, private communication), that 
is based on a larger and higher precision data set than Saar et 



al. J1998t ones and takes into account the color dependence, 
predict instead a jitter of about 2 and 4 m/s for components A 
and B respectively. 

To get a jitter of 8 m/s (required to explain the observed 
scatter), the chromospheric activity of HD 219542 B should be 
as high as log/?^^ ~ -4.75. Such a value seems to be excluded 
from the analysis presented in Sect. 12. 31 *. 

The observed radial velocity jitter of HD 219542 A is very 
low, at most 4 m/s, assuming that the 25% overestimate of 
errors of consecutive spectra holds also on longer time base- 



Note that our Rhk measurement was obtained in the season for 
which a larger radial velocity scatter is present. 
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line. On the basis of the results by Paulson et al. ilOOll for 
the Hyades, we do not expect that the colder component of a 
(coeval) pair shows larger jitter than the warmer one. However, 
the two stars may be on different phases of their activity cicle, 
so that the low jitter of HD 219542 A does not automatically 
guarantee the absence of moderate activity in its companion'^. 

Therefore, we can not exclude an activity origin for the ve- 
locity variations, taking into account that the line bisector, pho- 
tometry and equivalent width analysis do not have the precision 
required to evidence the variations causing a ~ 10 m/s radial 
velocity scatter 

7. Limits on planetary companions 

As shown in Sect. 14. 21 HD 219542 A does not show variations 
above errors and the statistical analysis did not reveal any sig- 
nificant periodicity in the velocities. The small velocity errors 
and the dense sampling of the observations allow us to derive 
firm upper limits for possible substellar companions. To this 
aim, we have developed a procedure that allow us to derive 
constraints for the masses of planets in circular as well as in 
eccentric orbit. 

The method is a Monte-Carlo extraction and is as follows: 
we set a dense grid of mass and semimajor axis values covering 
the ranges 0.05-4My and 0.01-3 AU. For each pair of mass and 
semimajor axis we randomly generated (with uniform distribu- 
tion) 10000 combinations of orbital phase To eccentricity e and 
longitude of the periastron co (the last two only for the eccentric 
case). We adopt the following ranges: e < 0.5625xlog(f )-0.12 
for f < 80 days (a conservative upper limit to the eccentricity 
of known exoplanets as a function of period) and < e < 0.95 
for longer periods < w < 360 deg; < Tq < Pmn.\- 

For each set of orbital parameters, we create a synthetic ra- 
dial velocity curve obtained by adding a Gaussian noise term 
(with the same rms as our RV-results) to a Keplerian orbit sam- 
pled at the epochs of our observations. 

We then applied the F test, and computed the probability 
that the synthetic and the observed data sets have the same vari- 
ance. If the probability is larger than 5% for more than 95% of 
the 10000 data sets for each mass-semimajor axis pair, than we 
took the planet with the corresponding mass and orbital semi- 
axis as compatible with the observed data. Otherwise, we con- 
sidered the planet as excluded by our data, since its presence 
would have created a detectable excess of variability. 

In this manner we obtain two regions in the mass- 
semimajor axis plane the first where planets can be excluded, 
and an area below our limits (see Fig. I15ll6t . where planets 
can still exist. It must be noted that the derived upper limits 
concern m sin /. 

Our method, based on the evaluation of the excess of ra- 
dial velocity variability caused by the presence of hypothetical 
planets, allows a complete exploration of the possible orbital 

If the planet around HD 219542 B does exist, the upper limit to 
the jitter estimated from residuals from the best fitting orbit, is 4.2 m/s, 
similar to HD 219542 A. 

This upper limit is justified on the basis on the detection of a planet 
with e=0.93 (HD80606b, Naef et al. i2001j . 
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Fig. 15. Limits on planets around HD 219542 A in circular 
(dashed line) and eccentric (continuous line) orbits derived us- 
ing the Monte Carlo procedure described in the text. 



parameters for eccentric orbits (the real case, since most of the 
known planets are in eccentric orbits). This is not feasible when 
the synthetic data are examined using methods like the peri- 
odogram analysis (see Endl et al.|2P02 1. 

Note that this procedure attempts to answers the question: 
would the presence of a planet more massive than this limit at 
the given separation have caused scatter incompatible with ob- 
served data? This is a somewhat different question from that 
answered by the usual procedure, that consider the detection of 
a planet with a high level of confidence (typically 99%) using 
procedures based on the periodogram analysis (see e.g. Endl et 
al.|2002| Cumming et al. 1999). These authors formulate the 
problem slightly differently: they determine the minimum am- 
plitudes for planetary signal at various separations, which can 
be detected by their data at a statistical significance of > 99%. 
The two approaches are rather different, because in principle, 
a planet could be present even if its detection is not significant 
at a high level of confidence. However, we repeated our proce- 
dure for the case of a Cen A, for which upper limits to planet 
masses have been obtained by Endl et al. ( 1200 It using the alter- 
native approach, and we find that our procedure (when assum- 
ing circular orbits) closely reproduces the same upper limits to 
planetary masses obtained in that paper. 

Fig.fBl shows the results for HD 219542 A. The Hmits for 
planets in eccentric orbit are less stringent than for planets in 
circular orbits, becuase of the flexibility of shape of the keple- 
rian curve for high eccentricties. Planets with msini > 2 Mj 
are excluded for separation smaller than 2.4 AU also 51 Peg- 
Uke planets with masses msini > 0.15 Mj are excluded. For 
orbital periods exceeding the baseline of our measurement the 
constraint becomes shallower. A peak corresponding to the 1 
year observing window can be noticed, especially in the eccen- 
tric case. Fig.ll6lshows the results for HD 219542 B. The limits 
are less stringent because of the larger scatter of the velocities. 
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Fig. 16. Limits on planets around HD 219542 B in circular 
(dashed line) and eccentric (continuous line) orbits derived us- 
ing the Monte Carlo procedure described in the text. 

8. Discussion 

The main results of this paper are the following: 

1. We detected a low-amplitude RV variation for 
HD 219542 B which could be due to a Saturn mass 
planet in close orbit; 

2. For component A, the star which shows signs of planetary 
accretion, we can place stringent upper limits on the exis- 
tence of orbiting giant planets. 

Taking into account the possibility that that the observed 
radial velocity variation of component B is caused by modest 
stellar activity, we think that a lengthy discussion on the im- 
plications of the presence of a planet around HD 219542 B is 
premature. 

However, we note that HD 219542 Bb, if it exists, should 
represent the third case of a planet in a binary system around 
the lithium poor component of the pair, after 16 Cyg B and 
HD 178911. In the case of HD 219542 and 16 Cyg"*, the star 
without detected planet shows also the signatures of accretion 
of planetary material. In these cases, the signatures of accretion 
might be anti-correlated with the presence of planets. 

HD 219542 B is by itself a metal-rich star, with metallicity 
typical for the stars with planets ([Fe/H]=H-0.20). If we take 
the abundances of HD 219542 B as the original values one 
for the system (i.e. assuming that no enrichment occuiTed for 
HD 219542 B) it implies that the system was originally very 
metal rich. This might suggest that accretion events are more 
likely on very metal rich stars, possibly because they form more 
planets and/or they have more massive disks. The full analysis 
of the SARG sample would tell us if such enrichment events 
are limited to very metal rich stars or they can happen also for 
metal poor stars. 

To our knowledge, a high precision differential chemical compo- 
sition analysis was not performed for HD 178911, but this would be 
complicated by the double line spectroscopic binarity of the primary. 



The lack of giant planets in close orbit around 
HD 219542 A leaves several possibilities for the origin of the 
observed signatures of accretion open. 

- Planet formation did not occur around HD 219542 A and 
a significant portion of the circumstellar disk fell into the 
stars possibly because of the dynamical perturbation of 
HD 219542 B. This possibility is neither excluded nor sup- 
ported by the possible stellar orbital parameters given in 
Sect. |3l If the planet around component B is confirmed, it 
would indicate that gravitational influence did not disrupt 
the planet formation process. 

- One or more giant planets orbiting around HD 219542 A 
were scattered toward the star by other planets in the system 
(Marzari & Weidenschilling 2002 ) or by the perturbation of 
HD 219542 B. In the first case, we would expect to find a 
planet in an external eccentric orbit. The continuation of 
our radial velocity monitoring could reveal that. 

- No giant planets formed in the system while terrestrial ones 
formed in the inner regions. Some of them may have been 
scattered into the star by mutual interactions. 

Our results exclude instead a scenario with a giant planet 
that formed near the ice boundary (4-5 AU) and migrated in- 
ward, causing the infall on the star of the inner part of the disk 
and/or terrestrial planets already formed (we would have de- 
tected such planet, unless the planet itself was engulfed by the 
star). 

The study of a few pairs clearly does not allow us to reach 
a firm conclusion about the frequency of pairs with chemical 
composition differences; a general picture is missing. How fre- 
quent are these pairs? Is there any link between the pollution 
phenomena and some properties like the physical separation of 
the components and the initial composition? How frequent are 
planets around the components of nearly equal mass binaries 
with separation 200-300 AU and what are their orbital charac- 
teristics? Are planets really missing around polluted stars? 

Our on-going survey including a sample of about 50 bina- 
ries might allow some deeper insight into these very exciting 
topics. 
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